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Abstract. We extend the covariance-matrix description of atom-light quantum 
interfaces, originally developed for real and effective spin-1/2 atoms, to include "spin 
alignment" degrees of freedom. This allows accurate modeling of optically-probed spin- 
1 ensembles in arbitrary magnetic fields. We also include technical noise terms that are 
very common in experimental situations. These include magnetic field noise, variable 
atom number and the effect of magnetic field inhomogeneities. We demonstrate the 
validity of our extended model by comparing numerical simulations to a free-induction 
decay (FID) measurement of polarized 87 Rb atoms in the / = 1 ground state. We 
qualitatively and quantitatively reproduce experimental results with all free parameters 
of the simulations fixed. The model can be easily extended to larger spin systems, and 
adapted to more complicated experimental situations. 



1. Introduction 

Atomic ensembles play an essential and growing role in quantum optics pQ, with 
applications in quantum networking [2], generation of optical quantum resources [HUE], 
and quantum-enhanced instruments [HJ El El IE] • The continuous- variable (CV) approach 
[9] efficiently describes experiments involving many quanta. The great majority of 
CV atomic ensemble experiments are performed with gaussian states, although non- 
gaussian atomic states JTU] are required for some tasks [11J. Gaussian states can be 
described very economically in terms of mean values and variances, whereas general 
states require a description exponential in the size of the system. Madsen and M0lmer 
[12] introduced a covariance matrix technique to describe collective spin variables in 
quantum interfaces between light and spin-1/2 ensembles. A wide variety of effects, 
including spatial and temporal inhomogeneities, loss, decoherence, atomic transport and 
projective measurements have been incorporated into this framework [13J. An important 
omission until now has been the description of larger-spin systems. Most atomic 
ensemble experiments are performed with alkali atoms and thus have ground-state spin 
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of at least 1, which implies tensorial light shifts [12J and generalized Faraday rotation 
effects [H]. An extension of the covariance matrix techniques to include these will 
allow statistical description of many large-spin applications, including quantum state 
characterization [15] and preparation [16], quantum chaos [17], optical magnetometry 
[SJ [T8J [HI] and quantum non-demolition measurement [201 121] • 

While several earlier works have applied the spin-1/2 framework to spin-1 or larger 
systems through the identification of a two-state "pseudo spin-1 / 2" sub-system, there are 
scenarios in which the dynamics naturally involves more than two levels, and requires a 
more expanded description. A clear example is a spin-1 or larger atom in the presence of 
both magnetic and optical fields. The magnetic field couples Zeeman states differing by 
Am = ±1, whereas the optical fields couple also states with Am = ±2 through tensorial 
light shifts. This system has been much studied using density-matrix approaches (15] . 
which describe fully the average single-atom properties but not the noise properties, 
which arise from correlations among the atoms. Consideration of the quantum noise 
in these systems motivates the current work, in which we extend the covariance matrix 
approach to spin-1 atoms. To the suite of techniques available for spin-1/2 ensembles 
[T2] IT5] , we add the ability to treat both vectorial and tensorial light shifts, technical 
noise due to uncertainty in the atom number, and dephasing due to magnetic field 
inhomogeneities. 

The paper is organized as follows: in section [2] we present the formalism, which 
employs eight orientation and alignment operators to describe the F = 1 collective 
atomic spin. In section [3] we analyze the spin dynamics in the presence of probing 
light and an external magnetic field, including coherent evolution, decoherence due to 
scattering of probe photons, and dephasing due to inhomogeneous magnetic fields. In 
section [4] we review the description of optical measurement within the covariance matrix 
formalism. In section [5] we describe the initial state including technical noise from 
uncertain atom number. In sections [6] and [7] we compare numerical results of our model 
with experimental data and identify the practical limits of the gaussian approximation 
in this system. 

2. Formalism 

We work with collective operators describing macroscopic numbers of particles, for 
which a continuous-variable description is appropriate. Throughout, we use the 
covariance matrix techniques [T2J [22], which are sufficient to describe the Gaussian 
states encountered in the great majority of continuous-variable experiments. 

2.1. Quantum polarization description 

Polarized light in continuous variables can be described with Stokes operators: 

S x = -aV^a S y = -aV^a S z = -a^a (1) 
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Where a = (a + ,a_) T and a + ,a_ are the annihilation operators for the left and right 
circular polarization and a x , cr y , a z the Pauli matrices. The S x , S y and S z Stokes 
operators represent respectively linearly polarized light horizontally or vertically, linearly 
polarized light on the ±45° direction, left and right circularly polarized light. They have 
the same commutation relations as angular momentum operators, [S x , S y ] = iS z and 
cyclic permutations. 



2.2. Description of spin-1 ensembles 



A single spin-1 atom is described by a density matrix with 8 degrees of freedom, which 
we express in terms of 8 single-particle operators \. These generalize the Pauli matrices, 
in the sense that they are traceless, Hermitian, and obey the orthonormality relation 
Tr(AiAj) = 25\y The first three operators are the components of the spin vector f, 
obeying [f x ,f y ] = if z . For illustration, we give the spin-1 matrix representation: 
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The others are rank-2 tensor operators, for which we use the symbol j, with components 
(again with the spin-1 representation for illustration): 
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Figure 1. Contribution dW/dXi of the / and j operators to the Wigner distribution 



representing the state p 



Radius indicates magnitude, warm 



(cold) colors indicate positive (negative) contributions. Axis markers indicate unity. 
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With quantization axis along z, j m describes the population imbalance between m F = 
and other states, while j X)V describe mp = ±1 coherences. represent = ±1 
coherences in other quantization axes. 

We note that the above operator definitions are spin-independent, and that the 
results in this manuscript follow from these operator definitions, not from the spin- 
1-specific matrix representations. As we shall see below, the most important coherent 
interactions: Larmor precession, Faraday rotation, and tensorial light shifts, can be fully 
described using the above operators, even for larger spin. The formalism developed here 
is thus applicable to some scenarios involving spin-3/2 and higher. Not all processes 
can be explained using just f and j operators, however. For example, with spin-2 
atoms modulated optical pumping in the presence of a B-field has been used to produce 
hexadecapole moment due to coherence between Zeeman states with Amp = 4 



2.3. Spin visualization 

From the orthogonality relation Tr[AiAj] 
p can be expressed as 



P 



-I + - AjAj 



2<5ij, an arbitrary single-atom density matrix 



(2) 
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Table 1. Commutation relationships for single-atom operators. 
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where Aj = TrfpAj. This suggests a visualization in terms of the spin Wigner distribution 
W(p), which is efficiently calculated as in Dowling, et al. [21]. In Figure [l] we show the 
differential contribution to the Wigner distribution dW/dXi. This shows, for example, 
that jxijy,jk and ji are related by spatial rotations while j m is distinct. 

2.4- Commutation relationships 

The operators f , j, and S have commutators given by 

[X a ,X b ] = icl aXb X k (3) 

where the A are f or j components and a sum is over k is implied. The structure 
constants c^ aXb are completely antisymmetric in the three indexes, and 



k 



c hh = 1 >J» = 2 J*3i = Jv3m = ^3 ( 4 n 

fz fz 1 Jm Ok v ' 

4f = 4f = 4f k = 4f k = !■ ( 5 ) 

To this we can add 

cf Sy = 1 (6) 

All structure constants not given above are zero. The commutators are given explicitly 
in Table [T] 

2.5. Collective spin operators 

To describe the ensemble we define collective operators. If AW describes atomic 
operators acting the z'th of N\ atoms, then A = A w . Explicitly for the vector (f) 
and tensor (j) collective spin operators: 

N A N A 

j = (7) 

1=1 1=1 

We note that these inherit their commutation relations from the microscopic operators: 
[A , Afe] = zc^ aAf, Afc. Finally, we define a phase-space vector to describe the state of the 
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-^pulses 

V=BffiF©j© S w (8) 

i=i 

where © indicates the direct sum and B is the magnetic field vector at the location 
of the ensemble. It should be noted that B is here a classical field, whereas the other 
components of V are operators. The inclusion of B in V allows for classical uncertainty 
about the field to be incorporated in a natural way into the calculations [25], as we 
describe below. 



2.6. Covariance matrix 

We work within the gaussian approximation, i.e., we assume that V is fully characterized 
by its average ( V ) and by its covariance matrix IV: 

r v = 1(V A V + (V A V) T ) - (V) A (V) (9) 

where A indicates the outer product. It will be convenient to define <5V = V — ( V ), 
the fluctuations of V about the mean, from which IV = |(<$V A 5~V + (<5V A <5V) T ) — 
(5V) A (5V). 



3. Dynamics 



We describe the most important dynamical effects for light-matter interfaces, namely 
the light-matter interaction that occurs when a pulse of probe light passes the ensemble, 
and the rotation due to a magnetic field. Both of these interactions produce coherent 
rotations, loss of coherence, and addition of noise. 



3.1. Light-atom interactions 

The light-atom interaction is described by an effective Hamiltonian which describes 
the dispersive effects of the electric dipole interaction in second order [2H1 EZl I2H]- In 
simulations and in practice, it is very convenient to employ a train of optical pulses for 
probing the ensemble. Even for continuous probing it is useful to treat the probe as a 
train of contiguous pulses, as this allows a course-grained description of the polarization 
evolution. During the time the m'th pulse is passing through the ensemble, the effective 
Hamiltonian is 

A(m) A(m) q(m) 

= G 1 ^F Z + G 2 (— — J x + — — J y ) (10) 

T T T 

where G\ and G 2 are coupling constants that depend on the geometry of the atomic 
ensemble and probe beam, the atomic structure, and the detuning from resonance [29J. 
We label the Stokes operators as S^ m '(t), m = 1 . . . A^ pu i ses where t m indicates the time of 
arrival at the ensemble. Thus S^ m \t < t m ) describes the polarization of the m'th pulse 
before entering the ensemble, while S^ m '(t > t m + r), where r is the pulse duration, 
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describes the polarization of the same pulse after exiting the ensemble. At time t m + r, 
the m'th pulse has left the ensemble (we assume the transit time is much less than r) 
and the change in the system is described by 

A(t m + r) = A(t m ) - zr[\(t m ), H™] (11) 

and 

S (m) (t m + r) = & m \t m ) -iT[& m Xt m ),Hg>]. (12) 

All other polarizations S^ n \ n ^ m are unchanged. These first-order difference equations 
will be accurate for sufficiently low-energy pulses, i.e., for small Gi j2 ( S ). For any 
given physical situation these conditions can be satisfied by subdividing long or high- 
energy pulses into sub-pulses with smaller r, at the cost of additional computation time. 
Similarly, if the full ensemble produces large rotations of S, the A can be subdivided as 
in [22]. In the simulations described below we subdivide S until the results converge. 

The evolution is compactly expressed in terms of a tensor H^ m ^ containing the 
coupling constants G\^ and Qf^q and the structure factors c^ Vj : 

* k 

V l (t m + r) = %{t m ) + VjCOi^VfcCO where (13) 
Hp = G 1 5f ) 4<?- + G 2 (sfcf^ + sfc^A (14) 

J K Vj V fc V Vj V k Vj Vfe J 

where 5^ 3 is 1 for Vj = Vj and for Vj ^ Vj, and we assume summation over repeated 
indices. 



3.2. Linearization 



The difference equations Eq. (13) are bilinear in the components of V. Although 
nonlinearity can in some cases lead to non-gaussian phase-space distributions [T71 [30] , 
in practice gaussian or near-gaussian distributions are far more common, and indeed 
producing measurably non-gaussian distributions is non-trivial [10J. This motivates a 
linearization of the above equations. Symbolically we write 

V(t m + r)= V(t m ) + V(t m ) ■ H< m > ■ V(0 (15) 

and separate V into the average V and the fluctuations <5V: 

V(t m + r) = V(t m + r) + SV(t m + r) 

= V(t m ) + V(t m ) ■ H< m > • V(t m ) + 5V(t m ) 

+ V(i m ) • H< m > • 5V(t m ) + 5V(t m ) ■ • V(t m ) 

+ 5V(t m )-H^-5V(t m ) (16) 

When the last term can be neglected, the dynamics reduce to 

V(t m + r) = V(t m ) + V(t m ) ■ H< m > ■ V(0 (17) 
= V(t m ) + V^-V(t m ) (18) 
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where U 1 -" 1 ^ = V(i m ) ■ U^ m \ which describes a nonlinear evolution of the average V, 
and 

5V(t m + r) = 5V(t m ) + V(U • H< m > • SV(t m ) 
+ 5V(t m )-n^-V(t m ) 
= T^-5V(t m ) (19) 

which describes a linear evolution of the fluctuations <5V in terms of the matrix 

Tt^St + VjHf^^V, (20) 

The covariance matrix T evolves as 

T(t m + r) = T^-T(t m )-[T^] T - (21) 
3.3. Optically-induced decoherence 

The above coherent rotations are necessarily accompanied by spontaneous scattering of 
photons, which produces decoherence of the atomic state [221 EZ] • An accurate treatment 
is possible [31], but requires a detailed accounting of the scattering channels. Here we 
describe a simple noise model, based on that of Madsen and M0lmer [12] . which agrees 
reasonably well with more sophisticated models [32] . 

We treat the decoherence as removal of a fraction 1 — X of the A? a atoms, followed 
by the re-addition of a fraction p of the removed atoms, with random polarizations. A 
common source of "loss" of atoms from the system is optical pumping into dark states 
[23 132] . In this article we take p — 1, i.e., no loss, which gives an upper bound on the 



introduced decoherence. As demonstrated in |Appendix A[ this alters the phase space 
distribution as 

A — > X\ (22) 

r A x 2 r A + x(i - *)A A r A + P {\ - *) jjiy A i 8 x8 (23) 

where A = (A), T A = 1(AaA + (AaA) t ) — (A)A(A) is the atomic part of the covariance 
matrix T v , I\ = |(A A A + (A A A) T ) — (A) A (A) is the single-atom covariance matrix, 
and the expectations are taken with respect to the average single-atom state p. T\ can 
be found from V as 



k 

M\f = l - (Tr[A fc {Ai, Xj}] - TrfAfcAijTrfAfcA,]) (24) 

where A = ( A ) is the single-atom average and { • , • } indicates the anti-commutator. 

The transformation of the atomic covariance matrix due to optically induced 
decoherence is then 

T A (t m + r) = Vi n) ■ T A (t m + r) ■ [V^Y + A^| m) (25) 
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where T^itm + r) is given by Eq. (21 ) and 

V { ™ ] = Xl 8x8 and (26) 
A^ m ) = X (i _ X)N A T X + p(l - X)^N A l 8x8 . (27) 

The fraction of atoms 1 — X that experience incoherent scattering of probe photons 
during a single pulse of duration r is given by X = exp(— 7/ 7 tol) , where r/ 7 is the 
probability for an atom to scatter a photon and n-^ is the number of photons in the 
pulse. A similar decoherence applies to the optical polarization [22J, where the 
fraction of photons scattered is 1 — e and e = exp(— tj^Na), leading to decoherence and 
noise terms 

^ m) = elsxs and (28) 
A^ m) ^ e (l-eHr s , (29) 

where Ts = (l/4)l 3x3 is the single photon covariance matrix. Because scattered photons 
are lost from the pulse, there is no photonic analogue of the last term of Eq. (27). In 
the scenarios described below the decoherence of the optical pulses is negligible, because 
the number of photons is much larger than the number of atoms, and we set e = 1. 
The full covariance matrix T then evolves as 

T(t m + r) = D< m > ■ ■ r(0 • [T< m )] T • [D^] T + N< m > (30) 

where = l 3x3 © © X>^ m) and N( m ) = 3x3 © N ( A m) © ■ 



3.4- Atom-field interaction 

In contrast to the atom-light interaction, the interaction of the atoms with the magnetic 
field is purely linear, so that precession by large angles can be described in a single step. 
At the same time, inhomogeneities in the magnetic field introduce dephasing, which 
requires a different description than given above for scattering-induced decoherence j|] 

To describe these effects, we first split the field into homogeneous and 
inhomogeneous parts as B(x) = Bo + B(x), with the assumption that |Bo| 3> |B(x)| 
for x within the atomic cloud. Specifically, if p a t(x) is the atomic density, we take 
B = f <i 3 xp at (x)B(x), the average over the atoms. 

Local rotation under B(x) can be described by the unitary operator 

U(pc,t)=exp[ijt\B(x)\f s (x)] (31) 

where /b(z) is the component of F parallel to the local field and 7 = iib9fI^ is the 
gyromagnetic ratio, where gp is the Lande factor and \Lb is the Bohr magneton. We 
note that f^ +2 (x) = f£(x) for n > 1, (true for any F — 1 spin component), so that a 
Taylor expansion of U gives 

C/(x,t) = 1 + /b(x) sm(~f\B(x)\t) + / B (x)[cos( 7 |i?(a;)|t) - 1] (32) 

I Here we assume that the magnetic fields are static, but time-varying fields can also be modeled, as 
described in Ref. [55] . 
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We can now see a qualitative difference between B||(x), the component of B(x) parallel 
to B , and Bj_(x) the perpendicular components. The effect of B||(x) is to change |B| 
and thus the precession frequency, causing a deviation from the behaviour under Bo 
that accumulates with time. In contrast, Bj_(x) principally changes the direction of /b, 
a non-accumulating effect. For this reason we focus on the effects of By(x). 

Using the Heisenberg equations of motion and the commutation relations of Eq. 
(|3|, we find dynamical equations 



d 
dt 



A(x) = -fi B g F \B(x)\A(x) ■ A(x) 
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A(x, t) = J2 e-^^Pi ■ A(x, 0) 



(34) 



(35) 



where w(x) = //s^|B(x)|, aj,Vj are the eigenvalues and eigenvectors of A with 
corresponding projectors Pi = v, A Vj. The eigenvectors are imaginary integers {a} = 
i{— 2, — 1, —1, 0, 0, 1, 1, 2}, indicating the half-period behaviour of some j variables. The 
macroscopic operators evolve as 



A(t)= I d 3 xp at (x)J2e~" { * )aitp i-Mx,0) 

= e ^° ait J d3x Pat(x)e-" (x)ait P 4 • A(x, 0) 



(36) 
(37) 



where w(x) = loq + u)(x). We now assume that A(x, 0) is uniform [§| so that the x- 
dependent terms can be collected as 7*j(t) = J Ap a t(x)e _IJ ' x '"*'. This gives 

A(t) = ^e-^^nMPi- A(0) 

i 

= V B (t)-T B (t)-A(0) (38) 

where Tsit) = Ylii e ~ u ° ait Pi describes the coherent evolution and T> B (t) = J2i r ii^)Pi 
describes the dephasing. 

In many situations it is reasonable to assume a Lorentzian distribution for Co. For 
example, in the highly-elongated trap described below, we observe an atomic density 



§ This assumption is clearly violated in some cases, for example in spin-echo experiments. However, it 
works remarkably well in many other situations. 
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p(z) well approximated by a Lorenzian p(z) = w/n(z 2 + w 2 ) where w ~ 48 mm is the 
full-width half-maximum (FWHM) extent of the ensemble. The decay term is then 



n(t) 



-wy\aidBu /dz\t 



or rj(t) = e l l Ti with 1/Tj = wy\aidBu/dz\, and 



(39) 



To preserve the uncertainty principle in the presence of the decay implied by Eq. 
(39), we must introduce a noise contribution N(t) to the covariance matrix [13]. The 
transformation due to evolution under an inhomogeneous field is thus 

r A (t) = v B {t) ■ r B (t) • r A (o) ■ j£{t) ■ vl(t) + j\f(t) (40) 

where the added noise obeys 

MB(t)>\iE'-iV B (t)-'E-Vl(t)\, (41) 
where I • I indicates the matrix absolute value and z'E*,- = 



[Aj,Aj]y and S' are the 
commutation matrices before and after the transformation |34|. Absent an indication 



that the dephasing should be extra noisy, we take the equality sign in Eq. (41). 



In addition to the dephasing of the atomic terms described above, the covariance 

ed 

\ 

7 ^(t) 7|B(x)|^ • 

V 



matrix T undergoes a coherent evolution as described by Eq. (21 ) in terms of the matrix 

/ 1,3X3 
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where the matrix A is described above in Eq. (34) and 
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/ 



Since the terms in TbW explicitly depend on the average atomic vector A(t) at 
time t, we update the covariance matrix in small time steps r < wo and keep track of 
the T(t) and V(£) at each step. The covariance matrix finally evolves according to 



T(t + r) = D B ■ T B ■ r(t) • [T B Y ■ [D B ] J + N B 



(44) 



where D B = 1 



3x3 



£>b © lsx3 and N 



3x3 



3x3- 



3.5. Combined effects 

To combine the coherent rotations with decoherence, we follow [13] and apply alternately 



Eqs. (17) & (21), Eqs. (22) & (30) and Eqs. (38) & (44) after subdividing the time 
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interval into sufficiently small steps r. Note that r is typically much smaller during 
the optical pulses than in the time between them, but in both cases we continue to 
subdivide the time evolution until the results converge. 

4. Measurement 

After passing through the ensemble, the optical pulses are typically measured by 
balanced polarimetry. This allows one component of S*™) to be detected. In 
the process, the pulse is absorbed by the detector and any possible information about 
the other components of is lost. We can describe the projective measurement by a 
unit vector p( m \ defined such that p^ m ^ • V = . 

The result of the measurement is random and gaussian-distributed, with a variance 
var(Sd c t) = p*™- 1 • T v •p (m * ) . The correlation of with other variables are given by 
the vector T v • p^. Finally, we can calculate the post-measurement uncertainty in V, 
in light of the measurement result and the known correlations, to find the posterior 
uncertainty 

(rypW)A(rypW) 
v v pw • r v • v {m) 
= r v - r v [n p(m) rvn p(m) ]( Mp )r^ (45) 

where the superscript ( Mp ) indicates the Moore-Penrose pseudo-inverse and IT p ( m ) = 
pM /\ pM i s a projector. 

It should be noted that there is nothing stochastic in the simulation method: only 
the average values, variances and covariances are computed. This again relies on the 
assumption that the fluctuations remain within a linear regime. 

5. Initial state and technical noise contributions 

A variety of initial states can be used with this approach, subject to a few limitations: 
The average value V must be physical, i.e., within the range of V for Na atoms and N L 
photons. The Robertson-Schrodinger uncertainty principle 

5A 2 5B 2 >\\([A,B\)\ 2 (46) 

places lower limits on the uncertainty implied by T v . Given any two scalar operators 
a = a-V, b = b-V 



1 



2 



5a 2 5b 2 = (a • T v • a)(b • T v • b) > - a t b jC y^v k . (47) 

In many situations it is appropriate to assume an initial state of the form R = p® NA 
for the atoms. For example, if optical pumping is used to initialize the state it is often 
reasonable to assume the atoms are independently pumped into a state p. This gives 
A = AATrfAp] and T A = NaT\. If p is a pure state, then R describes a coherent spin 
state (CSS). 
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We note an important subtlety about using CSSs for practical modeling: if p is 
a pure state, it is fully-polarized along some spin direction A p , with no uncertainty in 
that direction. This manifests as a zero eigenvalue of T\, and, if we naively apply the 
above, also in T A . This is unrealistic, however. In practice, N A is usually determined by 
a physical process with significant uncertainty, e.g. trap loading. Indeed, 5N A is often 
larger than the quantum variances N A Tx ~ A^a- This contributes a significant "atomic 
technical noise" to T A . As shown in Appendix B for a state of the form R = p® NA , 



where N A has average A^a and variance SN^_, 

T A = N A T X + 5N 2 A (\A\). (48) 

In most situations of interest the input light state is a coherent state. For example, 
in the experiments described below, we use an ^-polarized input, i.e. S = (Al/2,0,0), 
which has a covariance matrix Ts = diag(AT L /4, A^l/4, A^l/4). The input average 
magnetic field components B and covariance matrix Tb will depend on the modeled 
experiment. 

In the experiments described below the atoms and light are independently 
initialized, and we assume that the atom, light and magnetic field variables are initially 
uncorrelated. Thus the initial covariance matrix of eq. ^ can be written: 

r v = r B ©r A ©r s . (49) 

6. An example: Free-induction decay of collective atomic spin 

We illustrate the formalism described above with a simple example. We study, 
experimentally and theoretically, a paramagnetic Faraday rotation measurement of the 
free-induction decay (FID) of input F z and F y -polarized coherent spin states (CSS) 
precessing in a magnetic field. As described in Ref . [35] , neglecting tensorial light shifts 



described by the G2 term of Eq. 10, this leads to a measurable ploarization rotation 
angle 

U( t )\ = J ^ ByBz (1 " cos (7|5|t)e-* /T ) + B X \B\ sin( 7j Bt)e-*/ T ] (F s (0)> 

\ B \ 2 X I [B] + [B 2 X + B 2 y ] cos( 7 | J B|t)e-* /T ] (F z (0)) 

for the input F z and i^-polarized CSS respectively. The transverse relaxation time 
T = l/(wy\B'J) is due to the field-parallel gradient component Bt = d\B\/d z , and 
a Lorentzian distribution (full-width at half-maximum w) of atoms along z, the trap 



clXIS, clS described in |Appendix C[ Tensorial light shifts induce an additional nonlinear 
rotation of the atomic spins, as described in Refs [251 13T)] . This measurement can be 
used to estimate an unknown vector magnetic field, as in Ref. [35] , and is the basis of a 
proposal to prepare a planar-squeezed atomic spin state state [3"T] . 

Our experimental apparatus, illustrated in Fig. [2j has been described in detail 
elsewhere [8j |20j [38j [39] . In brief, we work with an ensemble of up to one million 
87 Rb atoms held in a weakly focused single beam optical dipole trap and probed by 
Us duration pulses of near-resonant light propagating along the trap axis and focused 
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Figure 2. (a) Experimental setup. PD: photodiode; L: lens; WP: wave plate; BS: 
beam splitter; PBS: polarizing beam splitter. Atoms are loaded into a single-beam 
optical dipole trap and probed with pulses of light propagating along the trap axis. An 
external magnetic field is applied to coherently rotate the atomic spins, and gradient 
field components dBi/dz are actively cancelled. The initial atomic state is prepared 
via optical pumping with circularly polarized light either propagating perpendicular 
to or along the trap axis, to prepare an i^-polarized or i^ 2 -polarized CSS, respectively. 
Also shown are experimental data of (b) the average signal <j> for an input F y - and 
i^-polarized CSS (blue and orange respectively), and (c) the evolution of the variance 
var(</>) for the same input states. Error bars represent ±lcr statistical errors. 



to a spot matching the radial width of the atomic cloud. This geometry produces a 
strong atom-light coupling, characterized by the effective on-resonance optical depth 
d Q = (o~o/A)Na, where o"o is the effective on-resonance scattering cross-section, A is 
an effective atom-light interaction area and is the number of atoms In this 

experiment we observe an effective optical depth do = 29.0(2). The optical readout can 
achieve projection- noise-limited sensitivity, calibrated against a thermal spin state [38J, 
and using dynamical decoupling techniques [20J has been used to demonstrate spin 
squeezing and entanglement-enhanced measurement sensitivity of the collective atomic 
spin [8]. 

For the FID measurement the atoms are initially polarized via optical pumping 
with a 5 |o.s duration pulse of circularly polarized light tuned to the / = 1 — > f = 1 
transition of the D 2 line. The optical pumping pulse propagates either along the trap 
axis to produce an F 2 -polarized CSS ((F x (0)) = or perpendicular to the trap to 
produce an F y -polarized CSS ((F y (0)) = Simultaneously with the optical pumping 
pulse, light tuned to the / = 2 — > f = 2 transition of the D 2 line is applied to prevent 
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atoms collecting in the / = 2 hyperfine state. This ensures that we prepare the initial 
CSS with good fidelity. We estimate an optical pumping efficiency of > 99% for both 
the .Fj-polarized and F y -polarized CSSs. The remaining unpolarized atoms make a 
negligible contribution to the observed var(0). 



full simulation 
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Figure 3. Illustration of the effect of technical noise terms and tensorial light shifts on 
the mean </> and variance var(</>) of the observed signal, (a) Tensorial light shifts rotate 
population out of F z and into the alignment variables J x y leading to a collapse and 
revival of the mean (j). The revival is accompanied by a 7r phase shift in the oscillations, 
as is evident in the simulation with G2 = 0, which removes the effect of tensorial light 
shifts, (b)-(d) The observed variance var(<^>) undergoes a more pronounced collapse and 
revival driven by the tensorial light shifts, which couple technical noise from the atomic 
variables A and the magnetic field covariance matrix Tb into the observed variable. 
The contribution to the observed signal in (b) is due to magnetic field noise. Technical 
noise in the atomic variables due to uncertainty SN^ in the atom number only adds 
significant noise only during the early stages of the evolution. This is illustrated more 
clearly in the magnified plots (c) and (d). The dot-dashed magenta line in plot (d) 
illustrates the quantum noise contributions due to the light (shot-noise) and atoms 
(projection-noise) . 



The atoms are then probed with a sequence of 1 u.s long pulses of linearly polarized 
light ((S x in) ) = ±N L /2) detuned about 1 GHz from resonance with the / = 1 —>•/' = 
transition on the D 2 line and propagating along the trap axis. After each pulse, 
the Faraday rotation signal = Sy^/Sx is recorded by a shot-noise limited 
polarimeter [40J . The pulses are either (a) all h-polarized and sent through the atoms at 
Pi = 10 \is intervals, or (b) pairs of alternately h- and v-polarized pulses separated by 
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3 us and sent through the atoms P2 = 20 us intervals. The entire experimental sequence 
is repeated up to 400 times to collect statistics. 

Typical experimental data are shown in Fig. [2j we plot the evolution of the average 
signal for an input F y - and F 2 -polarized CSS , and the evolution of the variance 
var(0) for the same input states. The observed signal initially oscillates at the Larmor 
precession frequency uq = 2n x 9.2 kHz with an exponential decay, as described by 



Eq. (50). On top of this there is a collapse and revival of the oscillations driven by 
tensorial light shifts, as described in Refs. [28l EE] , which rotate population out of F z 
and into the alignment variables J x>y at a rate uq 2 — {G2S X /2)(1/ Pi) = 2tt x 0.43 kHz. 
The revival is accompanied by a 7r phase shift in the oscillations (see Fig. [3] below for 
an illustration of this effect). The observed variance var(0) oscillates at a frequency 2coq 
and undergoes a similar but much more pronounced collapse and revival driven by the 
tensorial light shifts. Note that the observed variance in this experiment is dominated 
by technical noise; for our experimental parameters, the quantum noise contributions 
from both the atomic and light variables are <C 1 mrad 2 . The dominant technical 
noise contributions are due to uncertainty in the atom number 5N^, and magnetic field 
noise, described by the covariance matrix r#, coupled into the observed variable via 
tensorial light shifts. We illustrate the effect of these terms separately in Fig. [3j by 
running numerical simulations, as described in detail in Section [7| using experimental 
parameters from the data shown in Fig. [2] and setting variously the Tb, SN 2 and G2 
terms equal to zero. 

7. Simulation of free— induction decay 

In order to simulate the FID experiment, we need to estimate a number of experimental 
parameters, including the input state vector V and covariance matrix Tv, the coupling 



constants of Eq. (10), and the decoherence terms in Eqs. (26)-(29) and Eq. (39). The 
atom-light coupling constant G\ is calibrated in an auxiliary experiment as described 
in Ref . [39] . From this calibration we calculate the effective atom-light interaction area 



A, the coupling constant G 2 of Eq. ( 10 ) and the single-photon scattering probability rj 1 



used to describe the optically induced decoherence terms in Eqs. (26)-(29) (see Ref. [H] 
for details). 

We estimate the average magnetic field B for a single experimental data set by 



fitting the observed signal (4>(t)} with Eq. (50), as described detail in Ref. [35]. Since 



Eq. (50) neglects the effect of tensorial light shifts, we fit only the first 250 us of each 



measured signal to minimize the systematic error that this approximation introduces. 



From the fits we also determine the coherence time T that is used in Eq. (39) From 



the same fits we estimate the number of atoms (Fj(0)) = Na/2. Assuming an input 
state of the form R = p^ NA for the atoms, this specifies the inititial atomic state via 

Ji 

A = AaTt[Ap^J. The number of photons in each pulse is independently measured 
via a reference detector, as shown in Fig. [2] Together, these estimates specify the initial 
state of the vector V. 
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Statistics from the fits to Eq. (50) across the data set allow us to estimate the 
covariance matrix Tb- Similarly, we estimate the uncertainty 5N^ in the atom number 
from var(Fj) = vax(N A /2) (which includes contributions from variable trap loading 
and state preparation efficiency as well as measurement uncertainty). This allows us 
to estimate the initial covariance matrix T\ using Eq. (48). Since the measurement 
is shot-noise limited, the input light covariance matrix is Ts = N-^/A). 
Together, via Eq. (49), these specify the initial covariance matrix Tv(0). 

As an example, we give the experimental parameters in detail of the first example 
described in Section |6j shown in Fig. [2^b) and (c). For this experiment, we set the 
detuning of the probe to A = —700 MHz and probe with a sequence of 1 \xs long h- 
polarized pulses of light with on average = 7.2 x 10 6 photons per pulse at 10 |xs 
intervals. The measured average magnetic field was B = (11.98, —4.38, —4.01) mG, 
with a covariance matrix 

/ 0.202 0.0373 -0.048 \ 



V 



0.037 
-0.048 



0.201 
0.016 



0.016 
0.019 



mG' 



(51) 



7 



We estimate N A = 61.7 x 10 5 and AN a = 1.4 x 10 3 , giving an initial atomic vector 



for the F y -polarized input A(t 
matrix 



0) = (0, 1, 0, -0.5, 0, 0, 0, -0.29) x N A and covariance 
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(52) 



We further estimate a coherence time T = 360 ± 10 [is. The calibrated coupling 
constants were d = 1.7 ± 0.2 x 10~ 7 radians per atom and G2 = —7.5 ± 0.8 x 10 
radians per atom, and the atom-light scattering parameter i) 7 = 1.1 x 10 -9 . 

Together these parameters determine the initial state vector V(t = 0) and 
covariance matrix Tv(t = 0) as described in Section |5l the atom-light coupling constants 



of Eq. (10), the magnetic field components of Eq. (33), and the decoherence terms in 
Eqs. (26)-(29) and Eq. (39). With these input parameters, we then run numerical 
simulations of the evolution of the state vector V and covariance matrix Tv following 
the procedure described in Section [3j keeping track of V(t m ) and T v (t m ) at each time 
step. As in the experiment, we alternate between a time interval of 1 |js in which 
the light is present and an interval of 9 [\s with no light present. For the numerical 
calculations, we divide these intervals into 50 and 100 sub-steps respectively, which is 
sufficient to ensure numerical convergence of the results. The results of the simulations 
(dark blue curves) are plotted along with the experimental data (light blue circles) in 



Quantum atom-light interfaces in the gaussian description for spin-1 systems 



18 





60 




AH 


T3 
CO 


20 


DC 




E 











-20 




-40 








h (a) 

/l 


gauss 


experiment 






— simulation 








i i 







<|> experiment 
— simulation 



0.2 



0.4 0.6 
time (ms) 



0.8 



1.0 




0.4 0.6 
time (ms) 



1.0 



" (c) 


w 

CO . 

i 9 experiment 
simulation 







0.2 



0.4 0.6 
time (ms) 



0.8 



1.0 



1 (j> experiment 
simulation 




. (e) 


gauss 

— - 


experiment 


I A 




A A 


simulation 

ill Ui Ui 
A A /\j 








Iff IT TT 





140 
_ 120 
f 100 
1 80 

i 60 

40 
20 
0. 



0.2 



0.4 0.6 0.8 1.0 
time (ms) 



> 



(f) 


gauss 


1 


experiment . 




• »~ 

1 1 




- simulation - 




P 


Mi 













0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 
time (ms) time (ms) 

Figure 4. Comparison of experimental data (light blue circles) with numerical 
simulations (dark blue curves) of the mean c/> (a), (c) and (e) and variance var(0) 
(b), (d) and (f) of the optical rotation of the probing light. For (a) and (b) we set the 
detuning of the probe to A = —700 MHz and probed the atoms with a sequence of 1 ps 
long h-polarized pulses of light with on average 7V L = 7.2 x 10 6 photons per pulse at 
10 \is intervals. For (c) and (d) we set the detuning of the probe to A = —1.5 GHz and 
probed with the same measurement sequence. For (e) and (f) we set the detuning of 
the probe to A = —700 MHz and probed the atoms with a sequence pairs of 1 us long 
pulses with alternating h- and v-polarization, separated by 3 us and sent through the 
atoms 20 us intervals. Error bars represent ±ler statistical errors. See text for details. 



Fig. |4|a) and (b) for both the mean and variance var(0) with an input Fy-polarized 
atomic state. 

We observe excellent qualitative and quantitative agreement between the 
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simulations and the observed mean <f) and variance var(0) in the rotation angle with no 
further free parameters adjusted in the calculations. For the mean 0, we reproduce the 
observed behaviour in our simulations over the entire 1 ms of observation, as shown in 
Figs.[l|b) and (d). For the variance var(0) the quantitative agreement is initially good, 
but breaks down at longer times. This can be explained by the effect of the uncertainty 
in the magnetic field, described by the covariance matrix Tb, which eventually drives 
the atomic variables out of the gaussian approximation. 

Because magnetic precession is cyclic, not linear, an initially gaussian spin 
distribution will become non-gaussian due to uncertainty in the precession frequency. 
As some parts of the distribution precess faster than others, the distribution begins to 
"wrap around" the Bloch sphere, forcing a non-gaussian shape on the distribution. It 
is convenient to define a time r gauss = ir/^pABn), where AB» is the uncertainty in 
B along the average field direction, found using Tb of Eq. (51). r gauss indicates the 
moment at which the precession angle uncertainty becomes it. Perhaps surprisingly, 
the observed and predicted variances (see Fig|2|c) and (e)) agree very well up to 
t = r gauss w 0.45 ms, showing that the theory gives quantitatively accurate results 
even for significant departures from gaussianity. The agreement in the average values 
persists even for t > r gauss . This changes accumulates with time and becomes much more 
significant than the technical noise in the number of atoms so that while the atomic noise 
is important early on, the dephasing effect of magnetic field inhomogeneities becomes 
the dominant contribution for larger times. 

The collapse and revival of the oscillations in Figs. |4](a) and (b) is due to rotations 
driven by the tensorial light shift. The effect of the tensorial light shifts can be reduced 
either by probing further off resonance, or by probing the atoms alternately with h- and 
v-polarized pulses, as described in detail in Ref. [12] ■ We illustrate this in Fig. [4|c)-(d), 
where we compare these data to FID measurements (and simulations) made with two 
alternate probing strategies. In Fig. |4|c) and (d) we set the detuning to A = —1.5 GHz 
and repeat the single-polarization probing sequence. In Fig. |4|e) and (f) we set the 
detuning to A = —700 MHz and probe with pairs of alternately h- and v-polarized 
pulses separated by 3 \is are sent through the atoms 20 \is intervals (for clarity we plot 
only the h-polarized pulses). This results in the same total number of photons used per 
unit time as in the single polarization probing strategy. With these data we observe 
similar behavior in the mean and variance var(0) at both detunings, indicating the 
effective cancellation of tensorial light shifts in these experiments. We also observe 
non-zero minima in var(0) at all phases in the experimental data for t > 0.5 ms, which 
are not reproduced in our calculations. This may due to the finite sampling time of 
the measured data, which are not included within our calculations or in terms of the 
curvature of the Bloch sphere: for small uncertainty in precession angle A9 = t^pAB^, 
we will have AF Z oc A8dgF z . For flat parts of the curve F z {9) there will be zero d$F z 
and thus zero AF Z but for larger A8, we need to take into account higher derivatives, 
i.e. the curvature d^F z and higher. Finally we note that our model works always well 
for t < r gauss where we expect the Gaussian approximation to hold, as described above. 
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8. Conclusions 

We have presented a method for describing the quantum dynamics of spin-1 atomic 
ensembles, extending the method introduced for spin-1 /2 atomic ensembles by Madsen 
and M0lmer [12] and generalized by Koschorreck et al. [22J and Toth et al. [13]. Our 
approach, which explicitly includes the so-called "spin alignment" degrees of freedom, 
fits naturally with the light-matter interaction, which couples both spin alignment and 
spin orientation to the optical Stokes parameters. For spin-1 our description is complete 
within the gaussian approximation, while for larger spins it is still useful when octopole 
and higher spin moments can be neglected. We also include the important technical noise 
associated with magnetic fields and noise due to uncertain atom number, as typically 
arises due to stochastic trap-loading processes. Finally, we give explicit formulae for the 
noise introduced by spontaneous scattering during the optical probing process and due 
to dephasing in an inhomogeneous magnetic field. 

We have tested the model against experiment in a scenario involving all of these 
effects. We compute the evolution of the spin orientation average and variance for atomic 
ensembles with uncertain atom number, undergoing a combination of free-induction 
decay and alignment-to-orientation conversion in the presence of a noisy magnetic 
field. The simulation is compared to experimental observations made with a cold 87 Rb 
ensemble held in an optical dipole trap, probed by shot-noise-limited Faraday rotation 
with near-resonant light. We find good agreement within the gaussian regime. In 
addition to validating the model, the experiments provide a heuristic guide to the limits 
of the gaussian approximation in these systems. 

Given that most atomic ensemble experiments are performed with spin-1 or 
larger atoms, the technique described here will allow more accurate modeling of 
established quantum optical protocols, e.g. quantum memory [H], quantum non- 
demolition measurement [211 EH], dynamical decoupling [20], spin squeezing [5] and 
vector magnetometry [35J, as well as proposed applications such as generation of 
macroscopic singlet states [HH 05] and planar squeezed states [37] . 
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Appendix A. Proof of Eq. fl23| ) 

For brevity, we write the covariance as 

C(A, B) = Ua A B + (B A A)) - (A) A (B) (A.l) 

where the expectation is taken on the state of the ensemble. In this notation, 
T A ee C(A, A) where as above A = We assume that all atoms are statistically 

equivalent, so that (A« A\®) = (A (1) AA (2) ) for all % ^ j and (A«AA«) = (A^ AA«) 
for all i. We then have 

T A = $^C(A W ,A (J) ) = iV A C(A (1) , A (1) ) + N A (N A - 1)C(A (1) , A (2) ) 

(A.2) 

which we solve for C(AW, A^) to get 

Because of the symmetry, removing atoms does not change C(X^\ A^) or C(A^ 1 - ) , X^). 
If a fraction 1— X of the atoms is removed, the covariance matrix y^ Na ^ of the remaining 
atoms can be calculated as in Eq. ( A.2[ ), but summing i and j from 1 to XN A . We find 

t (xn a ) = xN A C(X {1) ,X m ) + XN A (XN A - 1)C(A (1) ,A (2) ) 
^ Na) X(XN a -1) fX(l-X)N* 

where 



N\ — 1 " V Na-1 



r x = C(A (1) ,A«) (A.4) 
is the single-atom covariance matrix. Dropping terms of order 1/N A and smaller, 

r (XN A ) = T (N A ) X 2 + TxX ^ _ ( A £j 

This accounts for the change in T A due to removing (1 — X)N A atoms. We must also 
account for the noise of returning a fraction p of these atoms to the ensemble in a 
decohered state. We assume they are completely random, and thus add the noise of a 
thermal state (i.e. variance /(/ + l)/3 per atom). 

r (XN A ) = t (Na) X 2 + VxX {i - X)N A + 2p{\ - X)N A l/3. (A.6) 
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Here 

k 

M^^^Tr[k{\M\ (A.7) 

where { ■ , • } indicates the anti-commutator. This can be shown using the expectation 
Tr[pAj] = Aj and the ort honor mality condition Tr[AjAj] = 25^ from which Tr[{Aj, Aj}] = 
4:Sij. We find the single-particle state p = |l + | J2i AjAj and the covariances 

cov(A i , Xj) = - ( {Xi, Xj} ) - ( Xi ) ( Xj ) 

= + hTr[X k {X h Xj}} - XiXj. (A.8) 

k 

Appendix B. Noise from uncertain atom number 

We consider the statistics of A for ensembles with N A atoms in a permutationally- 
invariant product state R( Na ^ = p 07VA , and taking a statistical average over Na- We 
indicate averages with subscripts, e.g. ( A ) R = Tr[A_R] indicates an expectation 
with respect to the state R, while ( A ) RN = ( Tt[AR( Na )] ) n indicates a statistical 
average of ( A ) R over the distribution of N A . Due to the structure of R( Na \ 
( )_r( jv a) = ( ^ } p = ( A ) p is independent of both k and of N A , so that 

N A 
k,l=l 

= N A ( XiXj ) p + N A (N A - 1) ( Xi ) p { Xj ) p (B.l) 

( AjAj ) r>Na = N A ( XiXj ) p + (5N 2 A + N 2 A - N A ) ( A t ) p ( Xj ) p 

(B.2) 

where 5N\ = ( N A ) — ( N A ) 2 indicates the variance. At the same time 

( Aj ) r>Na =N A (X l ) p (B.3) 

so that 

cov(A i , A ) r . Na = - ( AiAj + AjAi ) R Na - ( A, ) r>Na ( Aj ) r ^ Na (B.4) 

= iV A cov(A,, Xj) p + 5N 2 A (K) ( Xj ) p (B.5) 



and thus 



In terms of the single-atom covariance matrix T\ of Eq. (A.4), 

T = iV A r A + 5N 2 A {\ A A) (B.6) 
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Appendix C. Inhomogeneous magnetic fields 



The microscopic spin operators evolve as 



r(t) = R( Zi ,t)r(o) 

where P is the spin of the i'th atom with position Zi and 
R(z,t) = exp[ lF t\B(z)\A B ] 



(C.2) 



(C.l) 



where 



/ ~B Z By \ 

A B = B z — B x , 
V -B y B x j 



(C.3) 



is the generator of rotations about B and B 



B/|B|. Expanding the field as 



B(z) ~ B + (By + T5' ± )z + 0(z 2 ), where By is parallel to B and B' ± is perpendicular. 
We note that a change in the magnitude of B has an accumulating effect on the spin 
precession, i.e., the change in f grows with t. In contrast, a change in the direction 
of B has a fixed effect: From the perspective of the measurement, a rotation of B is 
equivalent to a rotation of both the initial state and the measured component F z . For 
small gradients <9^B <C B/Z atoms , where / a toms is the length of the cloud, we can ignore 
B' ± . This approximation, along with the fact that A n B +2 = —A B: allows us to write 



where ou(z) = 7f|B + zBy]. 

In our trap, we observe an atomic density p(z) well approximated by a Lorenzian 
p(z) = w/n(z 2 + w 2 ) where w is the full-width half- maximum (FWHM) extent of the 
ensemble. The collective spin F = J2 i f l then evolves as 



In the first term I + A 2 Bq describes a projector onto the direction of B . This is the 
steady-state polarization. The second line describes a decaying oscillation about B of 
the perpendicular components with a coherence time T = 1/(w7f|Bm|). 



R(z,t) ~ I + A Bo smtu(z)t + A 2 Bo [l - cosco(z)t] 



(C.4) 




v\A Bo sin w t - A 2 Bq cosa;ot)F(0) 



(C.5) 



(C.6) 



